Multifractal eigenstates of quantum chaos and the Thue-Morse sequence 
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We analyze certain eigenstates of the quantum baker's map and demonstrate, using the Walsh- 
Hadamard transform, the emergence of the ubiquitous Thue-Morse sequence, a simple sequence that 
is at the border between quasi-periodicity and chaos, and hence is a good paradigm for quantum 
chaotic states. We show a family of states that are also simply related to Thue-Morse sequence, and 
are strongly scarred by short periodic orbits and their homoclinic excursions. We give approximate 
expressions for these states and provide evidence that these and other generic states are multifractal. 

PACS numbers: 05.45.Mt, 05.45.Df 



The emergence of nonlinearity induced chaos from the 
hncar Schrodinger equation in the classical limit is a sub- 
tle and interesting phenomenon Q, |^ Q . The eigenstates 
of quantized chaotic systems are of current interest. They 
range from those that bear striking resemblance to classi- 
cal periodic orbits, and hence called "scarred states" 0, 
to those that are fruitfully described by the statistical 
properties of random matrices 0. Simple models, such 
as the logistic map, reveal the essence of classical chaos 
0. Comparable quantum models, such as the quantum 
cat map Q and the quantum baker's map |^, while pro- 
viding valuable insights, are less tractable. 

The classical baker's map Q, T, is the area pre- 
serving transformation of the unit square [0, 1) x [0, 1) 
onto itself, which takes a phase space point (<7,p) to 
where {q' = 2q, p' = p/2) if < g < 1/2 and 
[q' = 2(j - 1, p' = (p + l)/2) if 1/2 < q < 1. The stretch- 
ing along the horizontal q direction by a factor of two is 
compensated exactly by a compression in the vertical p 
direction. The repeated action of T on the square leaves 
the phase space mixed, this is well known to be a fully 
chaotic system that in a mathematically precise sense is 
as random as a coin toss. The area-preserving property 
makes this map a model of chaotic two-degree of free- 
dom Hamiltonian systems, and the Lyapunov exponent 
is log(2) per iteration. 

The baker's map was quantized by Balazs and Voros 
0, while Saraceno imposed anti-periodic boundary 
conditions. We consider exclusively the case N = 2^ for 
integer K. It has been known from early on that this 
case is special for the quantum baker's map 0, but this 
has not so far been made use of to analyze eigenstates 
more closely. The quantum baker's map, in the position 
representation, that we use in this Letter is: 
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where {GK)mn = {Pm\qn) = exp[-27ri(m + l/2)(n + 
l/2)/N]/y/N . The Hilbert space is finite dimensional, 
the dimensionality N being the scaled inverse Planck 
constant {N — 1/h), where we have used that the phase- 
space area is unity. The position and momentum states 
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FIG. 1: The participation ratio of the eigenstates of the quan- 
tum baker's map for A*" = 512 in the (a) position basis, and (b) 
Walsh-Hadamard basis. The states are arranged in increasing 
order of the latter number. 



are denoted as \qn) and \pm), where m,n — 0, ■ ■ ■ , N ~ 1 
and the transformation function between these bases is 
the finite Fourier transform Gk given above. The choice 
of anti-periodic boundary conditions fully preserves par- 
ity symmetry, here called R, which is such that R\qn} = 
\qN-7i-i)- Time-reversal symmetry is also present and 
implies in the context of the quantum baker's map that 
an overall phase can be chosen such that the momen- 
tum and position representations are complex conjugates: 
GKff' = 0*, if is an eigenstate in the position basis. B 
is an unitary matrix, whose repeated application is the 
quantum version of the full left-shift of classical chaos. 
There is a semiclassical trace formula, which, based on 
the unstable periodic orbits, approximates eigenvalues 
8]. 

Our first step in identifying special states is to Walsh- 
Hadamard (WH) transform Q eigenstates; if cj) is an 
eigenstate we study HK(j), where Hk = ®^ H, a in- 
fold tensor product of the Hadamard matrix H — 
((1, 1), (1, — 1))/-\/2. Very significant simphfications are 
seen in this basis for almost all states. In particular we 
analyze the number of principal components in the WH 
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FIG. 2: The intensity \{qn\4>tm)\'^ of the Thue-Morse state 
in the position representation. Shown are the states for A'^ = 
512 (plus) and A'' = 256 (circle) after the latter state was 
appropriately scaled. Also shown is the estimate for N = 
512 case (line). The left inset is the intensity of the WH 
tranform of the Thue-Morse state for A'^ = 512, while the 
right inset shows the same on a different scale. 



transform by means of the participation ratio PR// which 
is l/X^i \ {HK4')i\'^- In Fig. 1 we compare this participa- 
tion ratio with that calculated in the original position 
basis. For one remarkable class of states that is present 
for all iV, powers of two, PR// is the smallest and of the 
order of unity, for example when N = 1024, PR/f = 1.96 
for this state. This class of states seems to have been 
identified by Balazs and Voros Q, as the eigenvalue is 
very close to — i, however, the WH transform reveals 
their simplicity and beauty. In Fig. 2 are examples of 
these states, and their WH transform, where we also see 
their near self-similarity. We call these set of states, 4>tm 
the "Thue-Morse" states, as the principal peak in the 
Hadamard transform corresponds to the overlap with the 
final column (or row) of Hk which (apart from the fac- 
tor 1/V2^) is the K-ih generation of the Thue-Morse 
sequence. 

The Thue-Morse sequence is an example of a binary 
"automatic sequence" [l^, whose first few terms are 
{1, —1, —1, 1, — 1, 1, 1, — 1, • • • }. The finite generations are 
constructed as follows: start with Iq = {1} and generate 
ti = {1, —1} by appending to t^, where the overbar 
is multiplication by —1, and we continue to iterate the 
rule tk+i = {tk, tk\. At stage K we get Ik a string 
of length 2^, which we also treat as a column vector 
whose n-th element we denote as txin). The concate- 
nation rule above is then equivalent to the generating 
rules: txi'^'n) = tK-i{n) and txi'^n + 1) = —tK-iin). 
The Thue-Morse sequence considered as a word with two 
alphabets is cube free, that is no block repeats thrice con- 
secutively. This sequence occurs in numerous contexts 
|1C| . and is marginal between a quasiperiodic sequence 
and a chaotic one. The deterministic disorder of this 



sequence is relevant to models of quasicrystals [llj and 
mesoscopic disordered systems |12j , and as we show here 
quantum chaos. We believe that this is the first time it 
has emerged naturally in a quantum mechanical problem, 
rather than being assumed. 

That Ik is an approximate eigenvector of B follows on 
using the concatenation rule to get {BtK)i = 

^ N/2-l N/2-1 

where C; = (1 — i(— 1)'). For large the n sum 
sharply peaks at I — 2m, or 2m + 1 depending on if I 
is even or odd. Thus {BtK)2m ~ — i(2\/2/7r)tif _i(m) 
and {BtK)2m+i ~ «(2\/2/7r)ii<'_i(m), from which using 
the generating rules of the Thue-Morse sequence we get 
{BtK)m ~ — i(2V2/7r)tif (m). The eigenvalue we get 
has a modulus 0.9. We can alternatively calculate 
exactly Gk^k as a product if we use the fact that txin) 
is (—1)'^ where r is the number of times that 1 occurs in 
the binary representation of n. A short calculation gives 

{GKtK)m = 

K-l 

V2^(-l)™+i (i)^^+i W sin[7r2^-^(m-|-l/2)], 

a useful formula from which also we can see that Ik is 
an approximate eigenvector of B. 

For example when N = 2^ |(i9|0tm)|V^ ~ 0.74. Tak- 
ing into account symmetries helps us do much better. It 
is easy to see that RIk = (— l)^^^, and indeed the par- 
ity of the Thue-Morse state flips with each power of two. 
We note however that Gk tK 7^ t*K and therefore Ik 
as such does not have the correct symmetry properties. 
To take advantage of time-reversal symmetry we con- 
struct 00 = T^if + l*G^tKi which satisfies Gk4'o = 0o- 
It is straightforward to see that G^Ik is also an ap- 
proximate eigenstate of B. We determine numerically 
the complex constant 7 such that (j)Q best approximates 
the state (ptm- For example if A^ = 512, 7 is predom- 
inantly real and just 4>q — C (ig + Ggtg) is such that 
\{4>o\4>t7n)\'^ ~ 0.93, that is the state (j)tm is determined 
to more than 93% by the Thue-Morse sequence and its 
Fourier transform. C is a normalization factor and we 
have used that Gg^tg — Ggtg. Thus the determinisitic 
structural disorder of the Thue-Morse sequence is seen in 
a quantized model of a classically deterministic and fully 
chaotic system. 

We note that we can improve upon the simple ansatz 
00 above, by taking into account the second rung of peaks 
in the Hadamard transform of 0tm, K easily identifiable 
peaks (A^ = 2-^), as seen for example in the right inset of 
Fig. 2. These second rung of peaks also result from the 
Fourier transform of the Thue-Morse sequence, but not 

(r) 

exclusively. We introduce the notation t),, for an 2''-fold 
repetition of the Thue-Morse sequence of generation fc, 
for example column number 2^~^ — 1 of Hk is t^K-r- 
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The improved ansatz for the Thue-Morse state is then 

<t>L = 7o</'o + Ef=i(7fe + ll Gj})S^-^ t^U: where S is 
the shift operator, and we again determine the complex 
constants numcricaUy such that the ansatz is closest 
to the state. 

The shift operator S, acts on the position basis as 
S\qn) = \q2n) or \q2n-N+i) depending on if n < N/2 
or otherwise. We notice that S is "almost" B, only there 
is no momentum cut-off, as {pm\B\qn) = V2{pm\q2n) for 
n and m both < N/2 — 1. Infact S was proposed as 
the quantum baker's map by O. Penrose [13, however, 
apart from other issues, it does not respect time-reversal 
symmetry, in the sense that G^S*Gk 7^ S~^. If the po- 
sition state |g„) is denoted in terms of the binary expan- 
sion oi n — aK-iaK-2 ■ ■ ■ clq then S\aK^iaK-2 ■ ■ ■ oo) = 
\aK~2aK-3 ■ ■ ■o.Q'^K-i), hence = In. Such shift op- 
erators including momentum like bits were recently pro- 
posed as quantum baker's maps, however these also do 
not seem to respect time reversal symmetry |16|. We 
can find 7^ such that \{4>tm\4'trn)\^ ^ 0.998, this is shown 
in Fig. 2. The emergence of the shift operator is un- 
surprising as it commutes with the Hadamard transform 
Hk, and is close to the quantum baker's map B. The 
Thue-Morse sequence is an exact eigenstate of S and 
is hence a stand-alone state. 

It is known that the support of the spectral measure 
of the infinite Thue-Morse sequence is a multifractal , 
thus it comes as no surprise that (ptm has multifractal 
properties. We can study these either by comparing 
Thue-Morse states at different values of K or for a fixed 
K, by averaging over different sized bins. We perform a 
standard multifractal analysis 0| and see typical /(a) 
singularity spectrum. One simple diagnostic is the scal- 
ing of the inverse participation ratio (IPR). If 0„ denotes 
the position representation of the state 0, the IPR which 
is "^^Zq \<j)n\'^ scales as N~^^. When D2, the correla- 
tion dimension, is such that < D2 < 1, the function is a 
multifractal, between localized states {D2 — 0) and ran- 
dom states {D2 = 1). This IPR and multifractal analysis 
is carried out in the physically relevant position basis. 
Either of the two procedures gives for the Thue-Morse 
state the dimension D2 ~ 0.86, indicating its multifractal 
character. This scaling along with the full f{a) spectrum 
is shown in Fig. 3. 

The Thue-Morse state is by far not the only one influ- 
enced by the Thue-Morse sequence. We now construct 
a family of strongly scarred states consisting of K mem- 
bers. The PKh is low for these states also, for exam- 
ple when N = 512 PKh ^ 15. We identify this as a 
"family" based on the similarity of their PKh and the 
shared position of the K principal peaks in their WH 
transform. For up to = 2^'^ we have verified that this 
family exists and has K prominent peaks, at t^Li and 

all its — 1 shifts, 5''^t^l_i- Amongst this family, mem- 
bers are strongly scarred by period— 2, period— 4 and re- 
lated homoclinic orbits. Two examples of the members of 
this family along with their WH transforms are shown in 
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FIG. 3: (Color Online) The scaling of the IPR in the position 
basis as a function of A'^. Shown as points are the numer- 
ical data, while the smooth curves are the best fit straight 
lines. The cases are: (a) the Thue-Morse state, (b) a period— 2 
scarred state, (c) fixed-point scarred state after projecting out 
the uniform state, and (d) average over all the states. The cor- 
responding f{a) spectrums for the case A'^ — 8192 are shown 
in the inset. 
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FIG. 4: Two members of a family scarred by the period-2 
orbit for A'^ — 512. Shown are the intensities (in position 
basis) of the states in (a) and (c), while the corresponding 
WH transforms are shown in (b) and (d). 



Fig. 4. Taking time-reversal symmetry into account we 
may then write an approximate expression for this family 
of states based on K vectors, similar to the one above for 
the Thue-Morse state. For example for iV = 512 this pro- 
cedure enables us to reproduce the state strongly scarred 
by the period— 2 orbit up to more than 95% and most 
of the family to similar accuracy. The existence of such 
families may be indicative of systematics of periodic orbit 
scarring. 

This strongly scarred state and its family members also 
have multifractal properties, and indeed D2 ~ 0.8 for the 
period— 2 scarred state (Fig. 3). We have identified other 
families, but we now turn to the state that is strongly 
scarred by the fixed point (0, 0) and whose WH transform 
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FIG. 5: (Color Online) The density plots of \{qp\ip)\ on the 
unit phase-space square for the Thue-Morse state in the cases 
N = 128, 256, 512, and 1024 (left to right, top to bottom). 



is strongly peaked, indeed corresponding to t\^'^ /\/2^ ^ 
that is simply the uniform state. This also forms a series, 
for different N , whose eigenangles are close to zero. Since 
the Fourier transform of the uniform state given above 
tends to a constant (independent of N) at the origin (and 
at 1), we expect the scarred state to be dominated by this 
"Bragg peak" . In order to see any multifractal character 

even in this state we project tg^^ out of its WH transform 
and analyze the resultant state (after renormalizing) in 
the position basis. This shows that the "grass" even in 
this strongly scarred state is a multifractal with a value 
of D2 ~ 0.45, as shown in Fig. 3, where the /(a) spec- 
trum for the fixed point and period— 2 scarred states are 
also shown conforming to those expected of multifractals. 
For a gross measure of the multifractality of the entire 
spectrum, we averaged the IPR (in the position basis) of 



all the states and found a scaling with D2 ~ 0.9, that is 
shown in Fig. 3. 

Classical structures in quantum states are clearly seen 
in their phase-space or coherent state representations 
which is shown in Fig. 5 for the Thue-Morse states we 
have discussed above. Note the self-similarity of the 
four states shown here, the two different types of intri- 
cate structures arise from the alternating parity of these 
states. While dominated by the period-2 orbit and its as- 
sociated structures, we leave the details for a later work. 

We have merely skimmed the surface of the states in 
this Letter, highlighting the novel features found. The 
classical baker's map straightens out the intricate homo- 
clinic tangle of unstable and stable manifolds to lie along 
the q and p directions, which are also the basis in which 
we are viewing the quantum states. The multifractality 
and the simplicity of the states may originate in this, 
while in other models the complexity of the tangle may 
give rise to a "mixing" that makes them look random. 
That the dimensionality is a power of two is also a criti- 
cal requirement for the observations in this Letter. Qubit 
implementation of the quantum baker's map, such as the 
experimental work with NMR bits ^3 ^'^^ naturally of 
this kind. However self-similarity has been noticed indi- 
rectly earlier in the quantum baker's map in the corre- 
lations between the eigenstates and level-velocities 

III, 

even for N that are not powers of two. 
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